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The  purpose  of  this  paper  is  to  study  the  extendability  of  equilibrium  states  of  rodlike  nematic 
polymers  with  the  Maier-Saupe  intermolecular  potential.  We  formulate  equilibrium  states  as 
solutions  of  a  nonlinear  system  and  calculate  the  determinant  of  the  Jacobian  matrix  of  the 
nonlinear  system.  It  is  found  that  the  Jacobian  matrix  is  nonsingular  everywhere  except  at  two 
equilibrium  states.  These  two  special  equilibrium  states  correspond  to  two  points  in  the  phase 
diagram.  One  point  is  the  folding  point  where  the  stable  prolate  branch  folds  into  the  unstable 
prolate  branch;  the  other  point  is  the  intersection  point  of  the  nematic  branch  and  the  isotropic 
branch  where  the  unstable  prolate  state  becomes  the  unstable  oblate  state.  Our  result  establishes 
the  existence  and  uniqueness  of  equilibrium  states  in  the  presence  of  small  perturbations  away 
from  these  two  special  equilibrium  states. 
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Creative  Commons  Attribution  License,  which  permits  unrestricted  use,  distribution,  and 
reproduction  in  any  medium,  provided  the  original  work  is  properly  cited. 


1.  Introduction 

Liquid  crystal  polymers  (LCPs)  consisting  of  rigid  rodlike  macromolecules  in  a  viscous 
solvent  have  wide  technological  applications  [1-4] .  The  most  common  theoretical  framework 
for  modeling  nematic  polymers  is  to  represent  the  nematogenic  molecules  as  rigid  rods 
and  to  describe  the  ensemble  with  an  orientational  probability  density  function  (pdf)  [5]. 
Equilibrium  orientational  distribution  is  related  to  the  interaction  potential  by  the  Boltzmann 
relation.  Recently,  rigorous  mathematical  analysis  on  the  equilibrium  states  of  nematic 
polymers  has  garnered  serious  attention  from  mathematicians  led  by  Constantin  and  Titi 
[6-22] .  For  example,  various  proofs  for  the  axisymmetry  of  equilibrium  states  with  only  the 
Maier-Saupe  interaction  were  given  in  [8,  19]  for  the  2-D  case  and  in  [9,  18,  20]  for  the  3- 
D  case;  the  equilibrium  states  for  the  case  where  dipole-dipole  interaction  is  coupled  to  the 
Maier-Saupe  interaction  were  studied  in  [17, 21] .  Not  surprisingly,  many  mathematical  issues 
of  polymers  are  still  unexplored.  In  particular,  to  our  knowledge,  the  extendability  of  the 
equilibrium  states  has  never  been  attempted. 
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Abstract  and  Applied  Analysis 


In  our  previous  study  [21],  it  was  revealed  that  for  large  values  of  the  nematic 
strength  and  small  perturbations,  there  exists  at  least  one  solution  near  the  corresponding 
unperturbed  pure  nematic  solution.  In  [21],  we  concluded  the  existence  using  complicated 
free  energy  arguments.  Moreover,  in  [21]  only  existence  was  established  whereas  the 
uniqueness  was  not  addressed  at  all.  Our  goal  in  the  current  study  is  to  establish  rigorously 
the  extendability,  both  existence  and  uniqueness,  of  the  equilibrium  states  of  nematic 
polymers  in  the  presence  of  small  perturbations.  Notice  that  all  equilibrium  states  of  nematic 
polymers  are  axisymmetric  [9,  18,  20].  To  facilitate  the  discussion  of  extendability,  we 
introduce  two  new  variables  for  the  nonlinear  system  governing  the  equilibrium  states: 
one  variable  is  proportional  to  the  order  parameter  of  the  pdf;  the  other  variable  is 
proportional  to  the  biaxial  order  parameter,  which  measures  the  deviation  of  the  pdf  from 
being  axisymmetric.  One  advantage  of  this  approach  is  that  the  Jacobian  matrix  is  diagonal 
at  equilibrium  states.  We  will  show  that  the  Jacobian  determinant  of  the  nonlinear  system 
is  nonzero  everywhere  except  at  two  equilibrium  states.  These  two  equilibrium  states 
correspond  to  two  points  in  the  phase  diagram:  one  point  is  the  folding  point  where  the 
stable  prolate  branch  folds  into  the  unstable  prolate  branch;  the  other  point  is  the  intersection 
point  of  the  nematic  states  and  the  isotropic  states  where  the  unstable  prolate  state  becomes 
the  unstable  oblate  state.  The  extendability  of  the  equilibrium  states  except  these  two  states 
follows  immediately  from  the  implicit  function  theorem. 


2.  Nonlinear  system  for  equilibrium  states,  Jacobian  matrix,  and  its  determinant 

We  now  briefly  recall  the  mathematical  description  for  equilibrium  states  of  rigid  rodlike 
nematic  polymers.  The  orientation  direction  of  each  polymer  rod  is  denoted  by  a  unit  vector 
m.  In  this  study,  by  "pure  nematic  polymers,"  we  mean  the  case  where  the  Maier-Saupe 
interaction  potential  is  the  only  potential.  The  Maier-Saupe  potential  is  given  by 

Lt(m)  =  -b(mm)  :  mm,  (2.1) 


where  the  tensor  product  mm  and  the  tensor  double  contraction  A  :  B  are  defined  as 


mm 


m  in 

m2.m\ 

nijiri] 


m\m2  mirri3 

7772  7772  7772  7773  , 

77137712  77737773, 


A:Be  trace  (AB). 


(2.2) 


In  (2.1),  b  =  3N/2,N  is  the  normalized  polymer  concentration  describing  the  strength  of 
intermolecular  interactions,  and  (mm)  is  the  second  moment  of  the  orientation  distribution: 

(mm)  =  mmp(m)dm,  (2.3) 

J  ||m[|=l 


where  p(m)  is  the  orientational  probability  density  function  (pdf)  of  the  ensemble  of  polymer 
rods.  For  convenience,  potential  (2.1)  has  been  normalized  with  respect  to  kgT,  where  kg 
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is  the  Boltzmann  constant  and  T  the  absolute  temperature.  For  pure  nematic  polymers, 
equilibrium  states  are  described  by  the  Boltzmann  distribution  [5] : 


(2.4) 


where  Z  -  fs  exp[-LT(m)]dm  is  the  partition  function  and  S  is  the  unit  sphere. 

We  choose  the  coordinate  system  such  that  the  second  moment  (mm)  is  diagonal: 


'(ml)  0  0 

(mm)  =  0  (mf)  0 

.  0  0  (ml). 


(2.5) 


As  a  consequence,  the  Maier-Saupe  potential  can  be  written  as 


LT(m)  =  -b((m\)m\  +  (m^)m\  +  ( m\)m\ ). 


(2.6) 


The  most  significant  conclusion  for  pure  nematic  polymers  is  that  all  equilibrium  states  are 
axisymmetric  [9,  18,  20].  Since  not  all  equilibrium  states  of  a  perturbed  nematic  polymer 
ensemble  are  necessarily  axisymmetric,  to  study  the  extendability,  we  formulate  the  problem 
without  using  the  axisymmetry  so  that  non-axisymmetric  perturbations  are  allowed.  The 
axisymmetry  will  be  useful  in  our  analysis  because  the  extendability  is  determined  by  the 
Jacobian  determinant  evaluated  at  the  unperturbed  equilibrium  state,  which  is  axisymmetric. 
The  definition  and  derivation  of  the  Jacobian  matrix,  however,  require  non-axisymmetric 
formulation.  For  pure  nematic  polymers,  the  total  potential  is  completely  specified  by  the 
second  moment  (mm).  As  a  result  of  the  Boltzmann  relation  (2.4),  an  equilibrium  state  is 
completely  specified  by  the  second  moment  (mm).  Because  of  the  constraint  m\+m\+m\  -  1, 
an  equilibrium  state  is  completely  specified  by 


si  =  (ml),  s2  =  ( m\ ). 


(2.7) 


At  a  pure  nematic  equilibrium  state  without  perturbation,  we  have  Si  =  s2  or  Si  =  1  -  Si  -  s2 
(which  means  Si  =  S3)  or  S2  =  1  -  S]  -  s2  (which  means  s2  =  S3).  At  an  equilibrium  state,  we 
select  the  coordinate  system  such  that  the  equilibrium  satisfies  Si  =  s2.  In  terms  of  (si,  S2),  the 
Maier-Saupe  potential  can  be  expressed  as 


LT(m)  =  -b(s\m\  +  s2m\  +  (1  -  Si  -  s2)ml) 


(2.8) 


Consequently,  the  equilibrium  pdf  is  given  by 


(2.9) 


4 


Abstract  and  Applied  Analysis 


At  an  equilibrium  state,  (si,  S2)  satisfies  the  nonlinear  system: 

Si  -  (ml)  =0,  S2  -  (ml)  =  0. 


The  form  of  the  equilibrium  pdf  (2.9)  motivates  us  to  introduce  (771,772): 


m  =  b 


1  -  y  (si  +  S2) 


ri2  =  b 


(s2  -  Si) 
2 


(2.10) 


(2.11) 


Here  771  is  proportional  to  the  order  parameter  of  the  pdf  while  772  is  proportional  to  the  biaxial 
order  parameter,  which  measures  the  deviation  of  the  pdf  from  being  axisymmetric.  As  we 
will  see,  one  advantage  of  using  (771,  772)  is  that  the  Jacobian  matrix  is  diagonal  at  equilibrium 
states.  In  terms  of  (77 1, 772),  the  pdf  has  the  expression 


p(m, 


exp  [772(7712  -  777j)  +  7717773] 

Jg  exp  [772  (ml -  ml)  +  7717773]  dm 


(2.12) 


The  second  advantage  of  using  (771,772)  is  that  the  pdf  does  not  depend  on  b  explicitly.  Later 
on,  this  property  will  enable  us  to  write  b  as  a  function  of  r  which  is  the  equilibrium  value  of 
771.  At  an  equilibrium  state,  the  nonlinear  system  for  (771, 772)  is 

F1(rn,n2;b)  =  ^  \  (3(777^)  -  1)  =  0, 

b  1  (2.13) 

F2  (771, 772;  b)  =  j-  -  -  (m22  -ml)=0. 


In  this  paper,  we  study  the  extendability  of  equilibrium  states  of  nematic  polymers  in  the 
presence  of  small  perturbations.  We  consider  a  perturbed  version  of  system  (2.13): 


£1(771, 772;  b)  =  e\, 
£2(771, 772;  b)  =  e2. 


(2.14) 


Mathematically,  we  study  the  existence  and  uniqueness  of  solution  of  system  (2.14) 
near  a  solution  of  system  (2.13)  for  small  e\  and  e2.  According  to  the  implicit  function 
theorem,  the  existence  and  uniqueness  are  determined  by  the  Jacobian  determinant  of 
(£1(771, 772), £2(771, 772)).  To  calculate  the  Jacobian  matrix,  we  first  calculate  the  derivatives  of 
pdf  (2.12): 


0 

—p(m, 771,772)  =  (ml  -  (77i^))p(m, 771,772), 
— p(m,77i,  772)  =  (777^-777^- 


(m22  -ml))p(m, 771,772). 


(2.15) 


H.  Wang  and  H.  Zhou 


5 


Recall  that  at  a  pure  nematic  equilibrium  state  without  perturbation,  we  have  selected  our 
coordinate  system  such  that  q2  =  0.  Evaluating  the  partial  derivatives  of  F\  and  F2  at  the 
equilibrium  (q  1  =  r,  q2  =  0)  yields 


ni=r,n2=0 

*Zi=r,»Z2=0 


ni=r,ri2=0 

h=r,ri2=0 


^  -  <"*!)))  =  l-|var(m§), 

-^{ml(m22  -m\-  (m\  -  m\)))  =  0, 
-i((m2-m2)(m2-(ra2)))  =  0, 

\-\^m  2~  ml)  O2  -  m\  ~  ( ml  -  m\)))> 

22~m2i)2)- 


(2.16) 


Here  we  point  out  that  the  Jacobian  matrix  of  system  (2.13)  evaluated  at  an  equilibrium  state 
is  diagonal,  which  is  caused  by  the  axisymmetry  of  the  equilibrium  state. 

It  follows  that  the  determinant  of  Jacobian  matrix  of  system  (2.13)  at  the  equilibrium 
is 


/d(FuF2) 


H\=r,ri2=0 


13  2 

-  -  -  var  (m3 


(2.17) 


In  (2.17),  all  averages  are  evaluated  using  the  equilibrium  pdf 


p(  m,ni,n2)\m=rn^0 


—  exp  (rm\). 


(2.18) 


For  pure  nematic  polymers,  an  equilibrium  state  is  completely  specified  by  r.  The  governing 
equation  for  r  is  the  first  equation  of  (2.13):  Fi(r,0;b)  =  0,  whereas  the  second  equation  of 
(2.13)  is  satisfied  automatically  when  ^2  =  0.  We  introduce  a  new  function: 

HwV)  =  Fi(n,0;b)  =  |  -  1(3 (ml)  -  1).  (2.19) 


The  governing  equation  for  r  is  F(r;  b )  =  0.  It  follows  directly  from  the  definition  that  function 
F(ri ;  b)  is  related  to  the  (1, 1)  element  of  the  Jacobian  matrix  at  equilibrium  as 


d 

difr 


Fi(rii,ri2-,b) 


n\=r,tl2=0 


kin;b) 


n=r 


(2.20) 


As  we  will  find  later,  this  relation  is  a  key  tool  for  determining  the  sign  of  the  (1,1)  element  of 
the  Jacobian  matrix  at  equilibrium.  In  the  above,  we  have  used  variables  q  and  r  more  or  less 
interchangibly.  Later  on  when  necessary,  we  will  use  q  to  denote  the  independent  variable 
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of  functions  and  use  r  to  denote  the  particular  value  of  ri  that  satisfies  F(r;b)  =  0.  For  that 
purpose,  we  will  continue  using  these  two  variables. 


3.  Extendability  of  equilibria  of  nematic  polymers 

To  obtain  more  specific  properties  of  equilibrium  states  of  nematic  polymers,  we  rewrite  the 
governing  equation  for  r  using  spherical  coordinates.  We  select  the  7«3-axis  as  the  pole  of  the 
spherical  coordinate  system.  The  equilibrium  pdf  given  in  (2.18)  becomes 


p(<p,  0;  r) 


exp  (rcos2(/>) 

lyt  exp  ( rcos2(p )  sin  cp  dip ' 


(3.1) 


where  <p  is  the  polar  angle  and  0  the  azimuthal  angle.  Here  we  need  to  point  out  that  r  is 
not  the  radius  in  the  spherical  coordinate  system.  It  is  a  parameter  of  the  equilibrium  state. 
Substituting  equilibrium  pdf  (3.1)  into  (l/2)(3(?n2)  - 1),  using  a  change  of  variable  u  =  cos  (p, 
and  integrating  by  parts,  we  have 

(u2  -  u4)  exp  ( ru 2 
Jp  exp  ( ru2)du 


-(3(ot2}-1) 


Jo 


Thus,  F(ri;b)  defined  in  (2.19)  becomes  F(rp  b)  =  [1  /Z?  -  f(r0],  where  /(if)  is  defined  as 


fin)  = 


Jq  u2(1  -  u 2)  exp  ( riu2)du 

Jo  exP  ( W2)du 


(3.3) 


The  equation  for  r,  F(r;  b)  =  0,  becomes 


r 


1 

b 


fir) 


=  0. 


(3.4) 


Equation  (3.4)  has  a  trivial  solution  r  =  0  for  arbitrary  value  of  b,  which  corresponds  to  the 
isotropic  branch  of  the  nematic  polymer  phase  diagram.  At  r  —  0,  the  equilibrium  pdf  is 
uniform,  p(<p,  0;  0)  =  1  / 4jt.  It  follows  that 


W)  =  \. 

(m4)  =  (mf)  =  uAdu  = 

(mWj)  =  <mim3>  =  \  jQ  “2(!  -  u2)du  = 

var  (ml)  =  (m\)  -  ( m\  f  = 

(  O2  -  mlf)  =  (ml)  -  2  (mlm\)  +  (m\)  =  |  -  T  = 


1  f  h 


4 

15' 


(3.5) 
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ni=o,ni=o 


13  ,,.12 


1/1  =0,1/2  =0 


1  1  //  2  2\2\  1  2 
=  --»((»«2-m1)  )  =  ^-^5- 


b  2 


(3-6) 


Therefore,  when  bf  15/2,  the  Jacobian  matrix  is  nonsingular  and  thereby  the  isotropic 
equilibrium  is  extendable.  At  b  =  15/2,  the  isotropic  branch  intersects  with  the  nematic 
branch. 

The  remaining  solutions  of  (3.4),  if  any,  satisfy 

l~f(r)  =  0.  (3-7) 

In  [20],  it  has  been  shown  that  the  function  f(rf)  has  the  following  properties: 

(1)  /( 0)  =  2/15,  lim^+oo f(if)  =  0  and  lim ,,^_00/(^)  =  0; 

(2)  f(rf)  attains  its  maximum  at  i\  =  rf  =  2.1782879748  >  0,  where  the  maximum  is 
f  in*)  =  0.14855559992254  >  0; 

(3)  fin )  >  0  for  n<rf  and  f'(rf)  <  0  for  n  >  rf. 

A  graph  of  f(if)  is  shown  in  Figure  1(a).  Thus,  we  draw  conclusions  below  for  (3.7)  as 
follows. 

(i)  The  critical  value  of  b  is  b*  =  1/  f  (if)  =  6.7314863965. 

(ii)  For  b  <  b* ,  (3.7)  has  no  solution. 

(iii)  At  b  =  b*,  (3.7)  has  one  solution  r(b*)  =  rf  >  0. 

(iv)  For  b*  <b  <  2/15,  (3.7)  has  two  solutions:  ru(b)  >  rf  >  0  and  0  <  rM(b)  <  rf. 

(v)  At  b  =  2/15,  (3.7)  has  two  solutions:  ru( 2/15)  >  rf  >  0  and  r(2/15)  =  0. 

(vi)  For  b  >  2/15,  (3.7)  has  two  solutions:  ru(b)  >  rf  >  0  and  r/  (fa)  <  0. 

As  mentioned  before,  we  use  subscript  "U"  to  denote  to  the  "Upper"  part  of  the  phase 
diagram  where  r  >  rf,  subscript  "M"  to  denote  to  the  "Middle"  part  of  the  phase  diagram, 
where  0  <  r  <  rf,  and  subscript  "L"  to  denote  to  the  "Lower"  part  of  the  phase  diagram, 
where  r  <  0.  The  phase  diagram  for  nematic  polymers  is  shown  in  Figure  1(b).  In  the 
terminologies  of  nematic  polymers,  curve  segment  ru(b)  is  the  stable  part  of  the  prolate 
branch;  rM(b)  is  the  unstable  part  of  the  prolate  branch;  and  r/  (fa)  is  the  unstable  oblate 
branch. 

Recall  that  b  is  the  strength  of  the  Maier-Saupe  interaction,  which  is  proportional  to 
the  normalized  polymer  concentration  and  is  inversely  proportional  to  the  temperature,  r, 
the  solution  of  F(r;  b)  =  0,  has  the  expression  r  =  b-  (1/2) (3(m^)  - 1),  where  (1/2)(3(wt3)  - 1) 
is  the  order  parameter.  Thus,  from  physical  considerations,  it  is  desirable  to  use  b  as  the 
independent  variable  and  treat  r  as  a  function  of  b.  Flowever,  r(b)  is  a  multivalue  function 
of  b  and  for  b  <  b*  function  r(b)  is  not  even  defined.  Mathematically,  it  is  much  more 
convenient  if  we  use  r  as  the  independent  variable  and  treat  b  as  a  function  of  r.  b(r)  is  a 
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(a)  (b) 

Figure  1:  (a)  Graph  of  function  (b)  phase  diagram  of  nematic  polymers. 


single-value  function  of  r  and  is  defined  for  all  values  of  r  in  (-00,  +00).  Below,  we  will  adopt 
this  new  formulation  of  viewing  b  as  a  function  of  r.  In  the  new  formulation,  function  b(r) 
is  determined  from  (3.7)  as  b(r)  =  1  /  f(r).  In  terms  of  this  new  formulation,  the  stable  part 
of  the  prolate  branch  ru(b )  can  be  simply  represented  as  b(r)  for  r  >  if;  the  unstable  part 
of  the  prolate  branch  rM(b)  as  b(r)  for  0  <  r  <  if;  and  the  unstable  oblate  branch  ri  (b')  as 
b(r)  for  r  <  0.  Now  we  discuss  the  extendability  of  these  branches.  Using  relation  (2.20), 
F(rpb)  =  ij[l/b  -  f(rf)],  and  1/b  -  f(r)  =  0,  we  arrive  at 


_d_ 

drii 


Fi(ili/il2;b(r)) 


m=r,tii=0 


d_ 

dri 


F{ji;b(r)) 


v.~r 


1 

Hr) 


nf(n) 


n=r 


=  -rf(r) 


' 

>  0,  for  r  >  if 

<  <  0,  for  0  <  r  <  if 

>  0,  for  r  <  0. 


(3.8) 


To  study  the  (2,2)  element  of  Jacobian  matrix,  we  first  rewrite  {l/2)({m\-m7f)1)  by 
substituting  equilibrium  pdf  (3.1)  into  (1/2) -  mf)2)  and  using  substitution  u  =  cos (p, 


1 

2 


1  Jo  (1 


-  u2y  exp  ( ru2)du 
Jg  exp  ( ru2)du 


(3.9) 


Using  (2.16),  b(  1)  =  l//(r),  the  expression  of  /(r)  given  in  (3.3),  and  result  (3.9),  we  have 
d  F2(rii,ri2;b{r)) 

n\=r,ti2=0 


b(r)  2 

1  Jq  (1  —  u 2)  (5m2  -  1)  exp  ( ru2)du 
^  j^exp  (ru2)du 


(3.10) 


=  g(r). 
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It  is  straightforward  to  verify  that  g(0)  =  0.  Below,  we  want  to  show  that  g(r)  >  0  for  r  >  0 
and  g(r)  <  0  for  r  <  0.  To  facilitate  the  analysis  below,  we  write  g(r)  as  an  average 


(3.11) 


where  the  average  is  with  respect  to  the  pdf  p{u;r)  =  exp (ru2)/  J(j  exp (ru2)du. 

Lemma  3.1.  Function  g(r)  given  in  (3.11)  has  the  property  that  g(r)  =  0  implies  g'(r)  >  0. 

Proof  of  Lemma  3.1.  We  first  calculate  the  derivative  of  the  pdf:  ( d/dr)p(u;r )  =  (u2  - 
(u2))p(u;  r),  which  leads  to 


S'(r)  =  |<(!  -  m2)(5m2  -  1)(m2  -  <«2))) 

=  ^  ((1  -  “2)  (5u2  ~  if)  +  |(1  -  5(»2))g(r). 

Thus,  whenever  g{r)  =  0,  we  have  g'{r)  =  (1/20)  ((1  -  u2)(5u2  -  l)2)  >  0. 

Lemma  3.1  together  with  g(0)  =  0  leads  to  g(r)  >  0  for  r  >  0  and  g(r)  <  0  for  r  <  0. 
Combining  this  result  on  g(r)  with  result  (3.8)  on  -r/'(r),  we  arrive  at 

>  0,  for  r  >  if, 

<0,  for  0  <  r  <  if,  (3.13) 
<  0,  for  r  <  0. 


det 


9(Fi,F2) 


=  {-rf\r))g(r)=  < 


qi=r,rtz=0 


(3.12) 

□ 


Therefore,  we  conclude  that  all  nematic  equilibrium  states  (stable  or  unstable)  are  extendable 
except  for  the  equilibrium  state  at  (r  =  0,  b  =  15/2)  and  the  equilibrium  state  at  (r  =  if,  b  = 

b*). 


4.  Conclusions 

In  this  work,  we  studied  the  extendability  of  equilibrium  states  of  nematic  polymers  with 
the  Maier-Saupe  intermolecular  potential.  We  found  that  the  Jacobian  matrix  of  the  nonlinear 
system  is  nonsingular  except  at  two  special  equilibrium  states.  The  significance  of  this  result 
is  its  implication  on  the  existence  and  uniqueness  of  equilibrium  states  of  a  perturbed  system, 
in  the  neighborhood  of  the  unperturbed  equilibrium  states. 


Acknowledgments 

H.  Wang  was  partially  supported  by  the  National  Science  Foundation.  H.  Zhou  was 
supported  in  part  by  the  Air  Force  Office  of  Scientific  Research  under  Grant  no. 
F1ATA06313G003. 


10 


Abstract  and  Applied  Analysis 


References 

111  B.  Bird,  R.  C.  Armstrong,  and  O.  Hassaeer,  Dynamics  of  Polymeric  Liquids,  vol.  1,  lohn  Wiley  &  Sons, 
New  York,  NY,  USA,  1987. 

[2]  A.  M.  Donald,  A.  H.  Windle,  and  S.  Hanna,  Liquid  Crystalline  Polymers,  Cambridge  University  Press, 
Cambridge,  UK,  2nd  edition,  2006. 

[3]  S.  Hess  and  M.  Kroger,  "Regular  and  chaotic  orientational  and  rheological  behaviour  of  liquid 
crystals,"  Journal  of  Physics:  Condensed  Matter,  vol.  16,  no.  38,  pp.  S3835-S3859,  2004. 

[4]  A.  D.  Rey  and  M.  M.  Denn,  "Dynamical  phenomena  in  liquid-crystalline  materials,"  in  Annual  Review 
of  Fluid  Mechanics,  vol.  34  of  Annual  Review  of  Fluid  Mechanics,  pp.  233-266,  Annual  Reviews,  Palo 
Alto,  Calif,  USA,  2002. 

[5]  M.  Doi  and  S.  F.  Edwards,  The  Theory  of  Polymer  Dynamics,  Oxford  University  Press,  Oxford,  UK,  1986. 

[6]  P.  Constantin,  I.  G.  Kevrekidis,  and  E.  S.  Titi,  "Asymptotic  states  of  a  Smoluchowski  equation,"  Archive 
for  Rational  Mechanics  and  Analysis,  vol.  174,  no.  3,  pp.  365-384,  2004. 

[7]  P.  Constantin,  I.  Kevrekidis,  and  E.  S.  Titi,  "Remarks  on  a  Smoluchowski  equation,"  Discrete  and 
Continuous  Dynamical  Systems.  Series  A,  vol.  11,  no.  1,  pp.  101-112,  2004. 

[8]  P.  Constantin  and  J.  Vukadinovic,  "Note  on  the  number  of  steady  states  for  a  two-dimensional 
Smoluchowski  equation,"  Nonlinearity,  vol.  18,  no.  1,  pp.  441^43,  2005. 

[9]  I.  Fatkullin  and  V.  Slastikov,  "Critical  points  of  the  Onsager  functional  on  a  sphere,"  Nonlinearity,  vol. 
18,  no.  6,  pp.  2565-2580,  2005. 

[10]  M.  G.  Forest,  R.  Zhou,  and  Q.  Wang,  "Symmetries  of  the  Doi  kinetic  theory  for  nematic  polymers  of 
arbitrary  aspect  ratio:  at  rest  and  in  linear  flows,"  Physical  Review  E,  vol.  66,  no.  3,  Article  ID  031712,  9 
pages,  2002. 

[11]  M.  G.  Forest,  Q.  Wang,  and  R.  Zhou,  "The  flow-phase  diagram  of  Doi-Hess  theory  for  sheared  nematic 
polymers  II:  finite  shear  rates,"  Rheologica  Acta,  vol.  44,  no.  1,  pp.  80-93,  2004. 

[12]  M.  G.  Forest,  Q.  Wang,  and  R.  Zhou,  "The  weak  shear  kinetic  phase  diagram  for  nematic  polymers," 
Rheologica  Acta,  vol.  43,  no.  1,  pp.  17-37,  2004. 

[13]  M.  G.  Forest,  R.  Zhou,  and  Q.  Wang,  "Chaotic  boundaries  of  nematic  polymers  in  mixed  shear  and 
extensional  flows,"  Physical  Review  Letters,  vol.  93,  no.  8,  Article  ID  088301,  4  pages,  2004. 

[14]  R.  Zhou,  M.  G.  Forest,  and  Q.  Wang,  "Kinetic  structure  simulations  of  nematic  polymers  in  plane 
Couette  cells — I:  the  algorithm  and  benchmarks,"  Multiscale  Modeling  &  Simulation,  vol.  3,  no.  4,  pp. 
853-870,  2005. 

[15]  M.  G.  Forest,  R.  Zhou,  and  Q.  Wang,  "Scaling  behavior  of  kinetic  orientational  distributions  for  dilute 
nematic  polymers  in  weak  shear,"  Journal  ofNon-Newto?rian  Fluid  Mechanics,  vol.  116,  no.  2-3,  pp.  183- 
204,  2004. 

[16]  M.  G.  Forest  and  Q.  Wang,  "Monodomain  response  of  finite-aspect-ratio  macromolecules  in  shear  and 
related  linear  flows,"  Rheologica  Acta,  vol.  42,  no.  1-2,  pp.  20-46,  2003. 

[17]  G.  Ji,  Q.  Wang,  R  Zhang,  and  H.  Zhou,  "Study  of  phase  transition  in  homogeneous,  rigid  extended 
nematics  and  magnetic  suspensions  using  an  order-reduction  method,"  Physics  of  Fluid,  vol.  18,  no. 
12,  Article  ID  123103, 17  pages,  2006. 

[18]  H.  Liu,  H.  Zhang,  and  P.  Zhang,  "Axial  symmetry  and  classification  of  stationary  solutions  of 
Doi-Onsager  equation  on  the  sphere  with  Maier-Saupe  potential,"  Communications  in  Mathematical 
Sciences,  vol.  3,  no.  2,  pp.  201-218,  2005. 

[19]  C.  Luo,  H.  Zhang,  and  P.  Zhang,  "The  structure  of  equilibrium  solutions  of  the  one-dimensional  Doi 
equation,"  Nonlinearity,  vol.  18,  no.  1,  pp.  379-389,  2005. 

[20]  H.  Zhou,  H.  Wang,  M.  G.  Forest,  and  Q.  Wang,  "A  new  proof  on  axisymmetric  equilibria  of  a  three- 
dimensional  Smoluchowski  equation,"  Nonlinearity,  vol.  18,  no.  6,  pp.  2815-2825,  2005. 

[21]  H.  Zhou,  H.  Wang,  Q.  Wang,  and  M.  G.  Forest,  "Characterization  of  stable  kinetic  equilibria  of  rigid, 
dipolar  rod  ensembles  for  coupled  dipole-dipole  and  Maier-Saupe  potentials,"  Nonlinearity,  vol.  20, 
no.  2,  pp.  277-297,  2007. 

[22]  H.  Zhou,  H.  Wang,  and  Q.  Wang,  "Nonparallel  solutions  of  extended  nematic  polymers  under  an 
external  field,"  Discrete  and  Continuous  Dynamical  Systems.  Series  B,  vol.  7,  no.  4,  pp.  907-929,  2007. 


